Topological states in pyrochlore iridates: 
long-range anisotropy strongly competing with spin-orbit interaction 
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In the search for topological phases in correlated electron systems, iridium-based pyrochlores 
A2lr207 - materials with 5d transition-metal ions - provide fertile grounds. Several novel topologi- 
cal states have been predicted but the actual realization of such states is believed to critically depend 
on the strength of local potentials arising from distortions of IrOe-cages. We test this hypothesis by 
measuring with resonant x-ray scattering the electronic level splittings in the A= Y, Eu systems, 
which we show to agree very well with ab initio electronic structure calculations. We find, however, 
that not distortions of IrOs-octahedra are the primary source for quenching the spin-orbit interac- 
tion, but strong long-range lattice anisotropies, which inevitably break the local cubic symmetry 
and will thereby be decisive in determining the system's topological ground state. 



PACS numbers: 

Introduction - The discovery that in solids with weakly 
interacting electrons the spin-orbit interaction can in- 
duce insulating states of topological nature, prompts the 
intriguing question what kind of topological phases can 
develop in correlated electron systems, in which it is well- 
known that the Coulomb interactions between electrons 
can qualitatively change the character and properties of 
the electronic ground state [IHZ]- Intra and inter-orbital 
electron-electron interactions are very substantial in 3c? 
transition metal compounds such as copper-oxides and 
become progressively weaker when going to heavier tran- 
sition metal elements, i.e., dd and 5d systems, as the 
d orbitals become more and more extended. Yet the 
relativistic spin-orbit coupling (SOC), the root cause of 
a number of recently predicted topologically non-trivial 
electronic states, follows the opposite trend - it increases 
progressively. In 5d transition metal compounds, such as 
iridates, the interesting situation arises that these inter- 
actions meet on the same energy scale. The way in which 
the ensuing competition plays out delicately depends on 
the electronic hopping amplitudes, electron-electron in- 
teraction strengths and, in particular, the local energy- 
level splittings of the Ir bd electrons, which in principle 
allows the stabilization of novel electronic states of mat- 
ter such as strong or weak topological Mott states, an 
axion insulator or a Weyl semi-metallic state [IHS] . 

In the pyrochlore iridates of the type ^2^207 that we 
consider here, with A = Sm, Eu, Lu and Y, five electrons 
occupy the three Ir t2g orbitals, which reside on Ir^"*" 
ions inside corner linked IrOe octahedra, see Fig. [ij This 
leaves one hole in the shell to which thus 6 distinct 



quantum states (3 orbital and 2 spin) are available. 
When the local symmetry is cubic, so that it does not 
lift the degeneracy of the three t2g orbitals, the strong 
SOC splits the t2g states up into a pure j = 3/2 quadru- 
plet and a pure j = 1/2 doublet. The doublet is higher 
in energy and therefore accomodates the hole. Any ad- 
ditional crystal-field splitting, for instance of tetragonal 
or trigonal symmetry, lifts the degeneracy of t2g states 
and competes with the SOC, thus tending to quench the 
orbital moment. As the SOC is driving the formation 
of electronic states of topological nature, the outcome of 
this competition is decisive for the actual realization of 
any type of non-trivial topological ground state in py- 
rochlore iridates [SI [7] . 

Experimental RIXS results - We use Resonant Inelas- 
tic X-ray Scattering (RIXS) [8 to measure directly the 
energy of the different configurations of a hole in the 
Ir t2g shell of Y2lr207 and Eu2lr207 and to determine 
in that way the crystal-field energy splittings of these 
states. The single crystals of Eu2lr2 07 and powder sam- 
ples of Y2lr2 07 were grown by solid-state synthesis (see 
Supplementary Material). RIXS is a second-order scat- 
tering technique and can directly probe the electronic 
transitions within the Ir 5d manifold due to two sub- 
sequent electric dipole transitions (2p — ^ 5d followed by 
5c? — >■ 2p) [SI 12] . It is therefore a valuable technique for 
detecting transitions between crystal-field split Ir 5d lev- 
els and has been utilized for a variety of iridates [TUin^ . 
We determine the splittings by measuring the d-d tran- 
sition energies at the iridium L3 edge, with an incident 
energy, Ei = 11.217 keV, chosen to maximize the res- 
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FIG. 1: (a) Network of corner-linked Ir4 tetrahedra in 
A2lr207 pyrochlore iridates. (b) A central IrOg octahedron 
(Ir in dark blue, O in pink) and how it is connected to its six 
neighboring IrOe octahedra. The six trivalent yl-site atoms 
(green) form a hexagonal ring around the central octahedron 
and thus a manifestly non-cubic surrounding. 



onant enhancement of the spectral features of interest 
below 1.5 eV. It should be noted that varying the inci- 
dent energy did not result in a shift of any of the peaks 
but merely changed their intensities, a behavior associ- 
ated with valence excitations [8t il0nl2) . The experiments 
were carried out at the Advanced Photon Source using 
the 9ID beamline with a Si(844) channel-cut secondary 
monochromator and a horizontal scattering geometry. A 
spherical (1 m radius) diced Si (844) analyzer was used 
and an overall energy resolution of 300 meV (FWHM) 
was obtained. 

Due to experimental conditions, the spectra for 
Y2lr207 and Eu2lr207 were obtained at two different 
temperatures, 300 and 150 K. Since the thermal contrac- 
tion of Y2lr2 07 is extremely small, within tenths of a 
percentage between 300 and 150 K, and the local oxy- 
gen octahedra are unaffected by the temperature [13], 
we can conclude that this difference in temperature has 
minimum bearing on our results. In Eu2lr207 no no- 
ticeable momentum dependence was seen for these exci- 
tations, indicating their localized nature. For reference 
purposes, the Eu2lr207 spectra shown in Fig. 2 were col- 
lected at Q=(8.45,8.45,7). To quantitatively analyze the 
RIXS spectra, the various peaks were fitted with analyt- 
ical functions, as shown by the dashed lines in Fig. 2. 
The zero-loss peak in both spectra was fitted with a sin- 
gle Gaussian of FWHM ^ 300 meV, the instrumental 
resolution. The low-energy excitations, Ei and E2, had 
to be fitted with one Gaussian and one Lorentzian, re- 
spectively. The Lorentzian function was used in order to 
capture some of the tail on the high energy side (^ 1.5 
eV). The high-energy excitations (E^) were fitted with a 
Gaussian peak on top of a sloping background. Such a 
sloping background can come from charge transfer exci- 
tations which are expected to appear in this range. 
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FIG. 2: Resonant Inelastic X-ray Scattering spectra of 
Y2lr207 (top) and Eu2lr207 (bottom) at the iridium L-edge. 
Dashed curves are the result of fitting (see text). The E\, E2 
and E3 energy positions for ¥2^207 (Eu2lr207) are 0.53 eV 
(0.54 eV), 0.98 eV (1.00 eV) and 3.90 eV (3.70 eV). 

The RIXS spectra of both ¥2^207 and Eu2lr207 
(Fig. 2) show sharp features below 1.5 eV, correspond- 
ing to transitions within the Ir t2g levels, and a strong 
intense peak stretching from 2 to 5 eV that according to 
the calculations, see below, corresponds to d-d transitions 
between the Ir t2g and Cg levels. The position refinement 
for the three main peaks apart from the zero-loss peak 
results in =0.53 ± 0.05 (0.54 ±0.03), i;2 = 0.98 ± 0.05 
(1.00 ±0.03) and i;3 = 3.90 ±0.05 eV (3.70 ±0.05 eV) for 
Y2lr207 (Eu2lr207). Using a first, empirical ansatz we 
fit the energies of the three peaks to the eigenvalues of an 
effective single-ion Hamiltonian for the t2g orbitals of the 
form Hq = Xl ■ s + Al^ , where A is the SOC strength and 
A the t2g crystal-field splitting (see Supplementary Ma- 
terial) . The latter tends to quench the Ir orbital moment 
and is usually identified with distortions of the IrOg oc- 
tahedra [3 [11]. The fit of the RIXS data to the effective 
A-A model results in parameters A = 0.41 and A = 0.59 
eV for Y-227 and A = 0.43 and A = 0.54 eV for Eu- 
227. The value for A for both these materials agrees very 
well to values of 0.39-0.49 eV that can be extracted from 
electron spin resonance measurements on Ir''"'" impurities 
|15j . The magnitudes of A, 0.54-0.59 eV, however, are 
surprisingly large. To understand the size and elucidate 
the microscopic origin of this large crystal-field splitting 
-a crucial energy scale in determining the topological 
ground state of the electronic system- we have carried 
out a set of detailed ab initio calculations of the Ir d-level 
electronic structure on a series 5d'^ pyrochlore iridates: 
Sm-, Eu-, Lu- and Y-227. 

Ab initio calculation of d-d excitations - To investi- 
gate in detail the electronic structure and the essential 
interactions in the ^2^207 iridates, we rely on ab initio 
many-body techniques from wave-function-based quan- 
tum chemistry [16] . To calculate the local electronic mul- 
tiplet structure with this wave-function-based methodol- 
ogy we use embedded clusters of IrOg octahedra and a 
farther-out solid-state environment that reproduces the 
Madelung field in the cluster region, see Refs. [T71 fT5] . 

Multiconfiguration self-consistent-field (MCSCF) and 
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SmalrzOy 0.61 0.91 3.41-4.75 

EuzIrzOy 0.60 (0.54) 0.91 (1.00) 3.39-4.72 (3.70) 

LualraOr 0.57 0.92 3.49-4.88 

Y2lr207 0.58 (0.53) 0.94 (0.98) 3.48-4.84 (3.90) 



TABLE I: Relative energies of the split j = 3/2 states E\ and 
E2 as well as t^g to t2ge\ excitation energies in 227 iridates. 
The energies are the result of ab initio spin-orbit multiref- 
erence configuration-interaction calculations, see text. The 
experimental values for Eu2lr207 and Y2lr207 determined 
from their RIXS spectra in Fig. 2 are shown in bold, within 
brackets. 



multireference configuration-interaction (MRCI) calcu- 
lations [16] were carried out to this end on large clus- 
ters made of one reference IrOg octahedron, six nearest- 
neighbor (NN) ^-metal sites, and six NN IrOe octahedra. 
The solid-state surroundings were further modeled as a 
large array of point charges fitted to reproduce the crystal 
Madelung field in the cluster region. All calculations were 
performed with the MOLPRO quantum-chemical software 
[TO] . We used energy-consistent relativistic pseudopo- 
tentials for Ir and the A elements [^O^ES] and Gaussian- 
type valence basis functions. Basis sets of quadruple-zeta 
quality were applied for the valence shells of the central 
Ir^+ ions [inj and triple-zeta basis sets from the MOL- 
PRO library for the ligands [H] of the central octahe- 
dron and the NN Ir sites I^Hj- For the central Ir ions we 
also used two polarization / functions [20]. For farther 
O's around the NN Ir sites we applied minimal atomic- 
natural-orbital basis sets [25]. The / electrons of the 
Ln?'^ species were incorporated in the effective core po- 
tentials [211 [22 and the outer sp shells of the Ln?^ and 
ions were modeled with sets of [3s2p] functions [Hj- 
[23] , Crystallographic data as reported by Taira et al. 
[26] were employed. 

For the ground-state calculations, the orbitals within 
each finite cluster are variationally optimized at the MC- 
SCF level. All Ir functions are included in the active 
orbital space [IB], i.e., all possible electron occupations 
are allowed within the t2g set of orbitals. On-site t2g and 
t2g to Cg excitations are afterwards computed just for the 
central IrOe octahedron while the occupation of the NN 
Ir valence shells is held frozen as in the ground-state con- 
figuration. The MRCI treatment includes on top of the 
MCSCF wave functions single and double excitations [TB] 
from the O 2p orbitals at the central octahedron and the 
Ir 5(i orbitals. The spin-orbit interactions are accounted 
for in the computations as described in Ref. [?7] . 

For extracting the local Ir t2g splittings, the NN Ir^+ 
ions were explicitly included in a first set of MCSCF 
calculations. However, the presence of six open-shell Ir 
NN's makes the spin-orbit calculations cumbersome be- 



cause for seven 5(i^ sites (and nondegenerate orbitals), a 
given electron configuration implies 1 octet, 6 sextet, 14 
quartet, and 14 doublet spin states which further interact 
via spin-orbit coupling (SOC). To simplify the problem 
and reduce the computational effort, we therefore further 
replaced the six Ir'*+ NN's by closed-shell Pt''"'" ions 
and in this manner obtain the relative MRCI energies for 
the spin-orbit states presented in Tables I and II (labelled 
MRCI+SOC). 

The d-d level splittings calculated for cuprates [^BHSO] 
such as La2Cu04 and Sr2Cu03 and iridates [12 [H] such 
as Sr2lr04 by similar quantum-chemical techniques are 
in close agreement to the experimental values of these 
excitation energies. Also for the 227 iridium pyrochlores 
that we consider here, we observe that the calculated 
excitation energies of the Sd multiplets (the values of Ei , 
E2 and E^), see Table I, are in close agreement with the 
ones obtained from our RIXS experiments. 

Origin of d-d splittings - The very good agreement be- 
tween calculated and measured d-d level splittings forms 
the basis for a subsequent detailed analysis of the micro- 
scopic origin of the crystal-field splitting of the 5d levels. 
To this end, we first test the hypothesis that the split- 
ting A in the effective single ion A — A model is due to 
a distortion of the IrOe octahedra which lowers the local 
cubic symmetry to a trigonal (or even lower) symme- 
try. It turns out that the crystal structure of the 227's 
under consideration is fully defined by just three parame- 
ters: the space group number, the cubic lattice constant 
a and the fractional coordinate x of the O at the 48/ 
site [2S]. For a; = = 5/16, the oxygen cage around 
each Ir site forms an undistorted, regular octahedron. In 
our 227 Ir pyrochlores, however, x is always larger than 
Xc, which translates into a compressive trigonal distor- 
tion of the IrOe octahedra and hence a splitting of the 
5d electronic levels. To estimate how large the resulting 
trigonal crystal-field splitting is, we have performed a set 
of further ah initio calculations, but now for an idealized 
crystal structure with x = x^ and thus undistorted oc- 
tahedra. The resulting energies E2 in crystals with 
perfect octahedra, presented in Table II, are seen to be 
just slightly different from the excitation energies of trig- 
onally distorted systems. Thus local trigonal distortions 
of the IrOg octahedra are not the cause of the energy 
splitting A of the t2g orbitals. 

To understand the physical origin of the large splitting 
A one needs to consider the crystal structure of the A- 
227's in more detail. As shown in Fig. lb, the A ions 
closest to a given iridium ion form a hexagonal structure 
in a plane parallel to two of the faces of the IrOg octa- 
hedron. It is clear that the potential generated by these 
six A^^ ions breaks cubic symmetry even in the absence 
of trigonal distortions. The numerical results we obtain 
show that due to the large ionic charge at the A sites 
and the large spatial extent of the iridium 5d orbitals the 
anisotropic, planar Aq coordination actually has by far a 



4 



Undistorted Octahedron Without SOC 





E? 


El 




A" A 


Eu2lr207 
YalraOr 


0.67 
0.67 


0.89 
0.86 




0.34 0.27 
0.32 0.30 


TABLE II: 


Relative 


energies Ei 


and E2 of the split j = 



3/2 spin-orbit states for idealized crystal structures without 
trigonal distortion of the IrOe octahedra. A" and A are the 
splitting of the Ir t2g levels in quantum-chemical calculations 
without SOC in the undistorted idealized and experimental 
|26| crystal structures, respectively. 



pyrochlore materials [31]. Yet it is pertinent for many 
more crystal structures, in particular for layered quasi- 
2D perovskites or chain-like quasi-lD 5d transition metal 
systems P l fTH ] . 
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stronger effect on the t2g level splitting than the trigonal 
distortion of the IrOe cages. 

To bolster this conclusion, an underlying assumption 
- the interpretation of the excitation energies Ei_2 in 
terms of an effective A— A model - requires crosschecking. 
We can do so by calculating the t2g crystal-field split- 
tings directly by switching off the SOC in the quantum- 
chemical calculations. These calculations show that the 
ground-state wave function has aig symmetry. The 5d- 
level splittings without SOC's (see Table II) are only 
marginally different between the distorted, experimental 
crystal structure (A) and the idealized, undistorted one 
(A°). Besides thus confirming that the 5d splittings are 
due to anisotropy potentials beyond the nearest-neighbor 
ligand coordination shell, the values computed without 
SOC's also bring additional insight into the magnitude 
of these splittings. Whereas the effective values for A in 
the A- A model and RIXS data are 0.54-0.59 eV, the 
ab initio calculations without SOC yield A = 0.27—0.30 
eV, i.e., values about a factor two lower. This implies 
that in the pyrochlore iridates the effect of the relativistic 
spin-orbit coupling cannot be captured by oversimplified 
single- ion models like the A — A model. Additional de- 
grees of freedom, in particular the t2g — eg couplings, the 
presence of anisotropic spin-orbit interactions and 
the hybridization of the Ir t2g states with other orbitals 
in the system must be considered. 

Conclusions - Our experimental and theoretical 
demonstration of the presence of large Ir t2g splittings in 
pyrochlore iridates which arise from longer-range crys- 
tal anisotropy |3U [33] and that directly compete with 
spin-orbit interactions, does not only change the canon- 
ical view that only, potentially weak, local distortions of 
IrOg octahedra tend to quench the spin-orbit coupling in 
these iridium pyrochlores. The broader ramification is 
that the rather extended nature of the 5d wave functions 
renders the longer-range anisotropy fields to be of fun- 
damental importance throughout the hd transition metal 
series. Their physical effect is particularly striking when 
the local ligand field symmetry is high, e.g. cubic or 
close-to-cubic, but the point group symmetry in the crys- 
tal is lower. This includes, for example, osmium-based 
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SUPPLEMENTARY MATERIAL 



Crystal growth 



Mixtures of Y2O3 and Ir02 with purities of 99.99% 
were ground in stoichiometric molar ratios, pelletized, 
and then heated in air at 1000°C for 100 hours. The re- 
sulting material was reground, pressed into pellets, and 
resintered at the same temperature for an additional 150 
hours, with two intermediate regrindings. Powder X-ray 
diffraction measurements confirmed the phase purity of 
the resulting Y2lr2 07 sample to within the resolution 
of the measurement. A single crystal of Eu2lr207 was 
grown by the solid state synthesis method, as previously 
described in detail in Ref. [35] . A mixture of polycrys- 
talhne Eu2lr207 and KF (2N) were heated up to 1100 
°C and next cooled down to 850 °C at a rate of 2.5 
°C/h. Resistivity data on the resulting Eu2lr207 single- 
crystal shows almost metallic behavior, indicating that 
the sample displays slight excess of Ir, see the discussion 
in Ref. [Ml- 



Single- ion A — A model 



In the simplest approximation, the energies of the three 
main peaks in the RIXS spectra can be mapped onto 
the eigenvalues of a single ion Hamiltonian of the form 
Ho = Xl ■ s + All the basis of t2g orbitals, A is the SOC 
and A is a tetragonal or trigonal t2g splitting [6j [T4] . 
The eigenvalues of Hq and the splittings among the spin- 
orbit states are given by the following expressions: 
Eq = A(-l + 6- V9 + 26 + d^)/4:, El = A/2, and E^ = 
\{-l + 6+^/¥+W+¥)/A, where Eq is the energy of the 
ground-state i~l/2 doublet and 5 = 2 A/ A . If Ei~Eq and 
E2 — E0 are known from experiment, i.e., RIXS, simple 
estimates for A and A can be in principle derived from 
the relations above. In particular, A = 2(2£'io— £'2o)/(3— <5) 
and A = A(5/2, where S = -b-VW^, b=il+3a^)/il-a'^), 
a = E2o/{E2Q — 2Eio), E2o — E2~Eo, and Eio = Ei — Eo . 



